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ABSTRACT 

SCUSS is a M-band photometric survey covering about 4000 square degree of the South Galactic Cap, reach¬ 
ing depths of up to 23 mag. By extending around 1.5 mag deeper than SDSS single-epoch u data, SCUSS is 
able to probe much a larger volume of the outer halo, i.e. with SCUSS data blue horizontal branch (BHB) stars 
can trace the outer halo of the Milky Way as far as 100-150 kpc. Utilizing this advantage we combine SCUSS n 
band with SDSS DR9 gri photometric bands to identify BHB stars and explore halo substructures. We confirm 
the existence of the Pisces overdensity, which is a structure in the outer halo (at around 80 kpc) that was discov¬ 
ered using RR Lyrae stars. For the first time we are able to determine its spatial extent, finding that it appears 
to be part of a stream with a clear distance gradient. The stream, which is ~5 degrees wide and stretches along 
~25 degrees, consists of 20-30 BHBs with a total significance of around 6<x over the background. Assuming 
we have detected the entire stream and that the progenitor has fully disrupted, then the number of BHBs sug¬ 
gests the original system was similar to smaller classical or a larger ultra-faint dwarf galaxy. On the other hand, 
if the progenitor still exists, it can be hunted for by reconstructing its orbit from the distance gradient of the 
stream. This new picture of the Pisces overdensity sheds new light on the origin of this intriguing system. 

Subject headings: Galaxy: halo - Galaxy: structure - surveys 


1. INTRODUCTION 

The Pisces overdensity is a recent discovery, being one 
of the most distant substructures in the G alacti c halo. This 
substructure was first uncovered by ISesar et ah ( 2007 1 using 
RR Lyraes from the multi-epoch Sloan Digital Sky Survey 
(SDSS) Stripe 82 data, along with a number of other candidate 
overdensities. Located at RA ~354° and Dec ~0°, with a me¬ 
dian heliocentric distance of 81 kpc, this was named ‘Struc¬ 
ture J’. Subsequently, IWatkins et all (120091) used the same 
data, t his time analyzing the light-curves from lBramich et all 
( 2008 ). and independently found what appears to be the same 
structure, located at -20° <RA< 0°, -1.25° <Dec< 1.25° 
with a heliocentric distance of 80 kpc. They named it the 
‘Pisces overdensity’ after the constellation in which it is lo¬ 
cated. To confirm whether this photometric overdensity is 
truly a co herent structure, as opposed to a chance concentra¬ 
tion, iKollmeier et al . (12009*) obtained spectroscopy for eight 
RR Lyrae stars in the Pisces overdensity region and found 
five of them have a narrow range of velo cities, whic h sugg ests 
that the overdensity is genuine. Later lSesar et all (120Id) ob¬ 
served a further four RR Lyrae stars and, using the combined 
sample of 12 stars, confirmed the presence of a secondary ve- 
l ocity struc ture, which was tentatively found in the original 
IKollmeier et all (120091) sample. 

In terms of its extension, one of the most intriguing aspects 
of the Pisces overdensity is that it appears to be distributed 
over a large area on the sky. The original detection spans 10- 
15°, which at a distance of 80 kpc corresponds to a width of 
15-20 kpc. However, since the existing studies are based only 
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on the thin 2.5°-wide stripe Stripe 82 data, it is impossible 
to determine the full extent. This makes it hard to draw any 
firm conclusions as to the nature of the progenitor, beyond 
the fact that it appears more likely to be a tidally disrupted 
galaxy rather than an intact system. The full extent can only 
be determined by either deep wide-field photometric surveys 
or spectroscopy of faint halo tracers at distances consistent 
with the Pisces overdensity. 

In this study, we will use a deep photometric sky survey, 
called the South Galactic Cap u band Sky Survey (SCUSS), to 
investigate the extension of the Pisces overdensity. SCUSS is 
a M-band (wavelength ~3538A) photometric survey, covering 
around 4000 square degree (sq. deg.) in the South Galactic 
Cap region. 80% of the area overlaps with the southern SDSS 
data, but the SCUSS magnitude limit is 1-1.5 mag deeper 
than SDSS III DR9 M-band data (wavelength ~3551A). This 
deep SCUSS M-band data can probe a much larger volume of 
the halo, easily reaching beyond 100 kpc for Blue Horizontal 
Branch (BHB) stars. Since the Pisces overdensity is about 80 
kpc away, the depth and sky coverage of SCUSS is ideal to 
study the spatial extent of this structure. 

We organize the paper as follows: Section [2] introduces 
the SCUSS data; in Section [3] we select BHB stars from the 
SCUSS survey; in Section[4]we investigate the spatial extent 
of the Pisces overdensity using these BHBs; and in Section 
0we explore the distribution of giant branch stars; finally we 
present a summary and discussion in Section^ 

2. DATA OVERVIEW 

We mainly utilize the SCUSS M-bancQ and SDSS DR9 
single-epoch gri bands data in this study. SCUSS is an in¬ 
ternational collaboration between the National Astronomical 
Observatories of China (Chinese Academy of Sciences) and 
Steward Observatory (University of Arizona, USA). This sur¬ 
vey has imaged -4000 sq. deg. of the Southern Galactic 

6 In the following, unless otherwise stated, when we discuss the M-band 
magnitude this corresponds to the SCUSS magnitude. 
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Cap, with galactic latitude b < -30° and equatorial latitude 
Deo -10° at a wavelength of u ~3538A. The survey was 
carried out using the 90Prime imager of the 2.3m Bok tele¬ 
scope at Steward Observatory on Kitt Peak between 2010 and 
2013. The detector consists of an array of four 4Kx4K (64- 
megapixel) CCDs and the field of view is 1.08° X1.03°, result¬ 
ing in a pixel resolution of 0.454". The typical seeing during 
the observing period is about 2.0" and the exposure time of 
each image is 300s. To obtain the required magnitude limit 
two dithered exposures for each field were taken and if any 
fields did not meet the quality requirement (for example due 
to bad weather), extra exposures were added. Mo re de tails 
abo ut the SCUSS su rvey can be found in IZhou et alJ (l2015h 
and lZou et al.1 (l2015h . 

The SCUSS observing strategy yields a significant increase 
in depth compared to single-epoch SDSS data. We compare 
these two surveys in Figure [I] where SCUSS co-added PSF 
magnitudes are compared with SDSS PSF magnitudes. In this 
figure we divide the whole SCUSS footprint into 1.08° x 1.03° 
bins and analyze the magnitude distribution (i.e. luminos¬ 
ity function) of point sources in each bin, where star-galaxy 
separation is based on the SDSS classification. The lumi¬ 
nosity function turns over at faint magnitudes as the survey 
reaches its detection limit and becomes incomplete. In our 
case the limiting magnitude is defined as the maximum mag¬ 
nitude of a star with a given error. For both SDSS & SCUSS 
datasets we consider detections with (/-band photometric er¬ 
ror less than 0.2 mag (5cr). This procedure is illustrated in 
the top panel of Figure [Q For this example field (located at 
RA = -41.5°, Dec = -8.5°) we find that the SCUSS magni¬ 
tude limit is around 23.3, which is around 1 mag deeper than 
SDSS. The bottom panel of this figure shows the distribution 
of magnitude limit across the SCUSS footprint. In general 
the SCUSS fields are 1 to 1.5 mag deeper. As our BHB clas¬ 
sification requires ugri photometry, we have only analyzed 
the 3400 sq. deg. of the SCUSS footprint that overlaps with 
SDSS. As can be seen from this figure, the overlap region 
is not contiguous due to the incomplete SDSS coverage for 
RA > 30°. The resulting depth varies somewhat across the 
footprint, but for most fields the limiting magnitude lies in the 
range 23 to 24 mag. Completeness for a given field is related 
to the turn-over in this luminosity function, which is typically 
0.5 mag brighter than the limiting magnitude. We have also 
measured the location of this turn-over for each field and find 
that only 2 per cent have a value brighter than 21.5 mag. For a 
BHB star this magnitude corresponds to a distance of around 
100 kpc, so we are confident that our BHB samples should be 
reasonably complete to 100 kpc and are able to probe signifi¬ 
cantly further than this (albeit at lower completeness). 

3. BHB SELECTION 

The class of A-type stars, which includes BHBs, can be 
selected using a color-color box in the space of (u - g) o vs 
(g - r) o. In the following work all magnitudes are corrected 
for extinction ( labe led with a subscript ‘o’) using the maps of 
ISchlegel et al.l(ll998l ) and ad opting the reddening conversions 
from lSchlaflv & FinkbeineJ (1201 lb . For SCUSS objects with 
16 < «o < 22.5 and u gri error less than 0.2, the A-type sta rs 
are selected following (lYannv et al.ll2000tiSirko et~aill2004l) . 

0.9 < (m - g) o < 1.4, - 0.3 <(g- r) 0 < 0.0. (1) 

The color-color criterion is efficient at removing contamina¬ 
tion of non-A-type main-sequence stars, white dwarfs and 




Fig. 1. — An estimation of the SCUSS magnitude limit. The top panel shows 
an example of how we define the magnitude limit. The red histogram shows 
the luminosity function in one SCUSS field (located at RA = -41.5° and Dec 
= -8.5°), while the blue histogram is for SDSS data in the same field. We use 
PSF magnitudes and only consider point sources with w-band error less than 
0.2 mag (5cr). We take the maximum magnitude (vertical dashed lines), as 
our definition of the magnitude limit. The bottom panel shows the resulting 
SCUSS magnitude limit for each field in the 3400 sq. deg. overlap region of 
SCUSS and SDSS. 

most of the quasars. To further eliminate quasars from the 
A-ty pe samp le we use the following color cut (similar to 
IDeason et al.ll2012l) . 

Cg ~ r )o > 0.6164(g - i) 0 - 0.016. (2) 

This cut removes many spurious points at the left and right 
edges of the (u - g)o vs (g - r)o diagram, making the BHB 
claw look more distinct. By cross-matching with SDSS 
spectroscopy we found that many of the points removed by 
the color cuts in Equation Q are indeed quasars. In order 
to check how many potential quasars are left in the result¬ 
ing sample, we apply the XDQSCQ technique (IBow et all 
1201 lalfbh . which has been designed to use for the SDSS quasar 
targeting selection. We found that the number of quasars 
(with probability greater than 0.8) left in our A-type sample 
is negligible, with a contamination fraction of less than 2%. 
IDeason et al.l (120121) also demonstrated the efficacy of their 
gri cut, finding that their spectroscopic sample of 19 faint 
(20 < g < 21.5) BHBs candidates contained no quasar con¬ 
tamination when this cut was applied. It is possible that there 
could be some variable star contamination (e.g. RR Lyraes) 
in the current sample, although it has been estimated that the 
contamination should only be around 5% and hence this will 

7 http://www.sdss3.org/svn/repo/xdqso/tags/v0-6/doc/build/html/index.html 
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have little bearing on our results dDeason et al.ll2Q lTh . Also 
note that since uq < 22.5, which is equivalent to a depth of 
r as 21.5 , the SDSS star-galaxy separation should be robust 
dLupton et al.ll200il see also the SDSS websit^l). 

To select BHB stars from A-type sample, we need to sep¬ 
arate them from Blue Stragglers (BS). BS stars also lie in 
the A-type regime and so are included within our (u - g)o 
vs (g - r )o selection (Equation Q}. To discriminate between 
BHB/BS stars, the ideal approach is to analyze the stars’ spec¬ 
tra, because the two gro ups have different surface gravities 


and spectral line profiles (iKinman et al. 

1994j (Clewley et al. 

2002; Sirko et al][2004; Xue et al. 2008[ 

Deason et al.l 2011[). 


BHB and BS stars and pinpointed the loci of the two popula¬ 
tions in the (n - g)o vs (g - r) o color plane. From their color- 
color diagram, the two groups have a distinct_division along 
the (u - g) o direction (see Figure 2 of iDeason et al.ll2011l) . 
This phenomenon is a result of the different surface gravi¬ 
ties of BHB and BS stars; since the n-band filter is located 
blue-ward of the Balmer discontinuity, the (m - g) o color char¬ 
acterizes the strength of the Balmer jump , a q uantity which 
is sensitive to the surface gravity. IDeason et al . (1201 ll) found 
that the (u - g)o distributions of BHB/BS stars can be fit us¬ 
ing one Gaussian for each population, with the center of each 
Gaussian varying with (g - r) () . By dividing the color-color 
plane into several slices in (g - r) o, and modeling the BHB/BS 
distributions of each slice with a two-Gaussian function, one 
can easily obtain a boundary line for BHB/BS stars. 

We classify BHB/BS stars using the above method. 
However, we cannot directly apply the boundary line of 
IDeason et al.l (1201 ll) as the u filters of SCUSS and SDSS have 
slightly different profiles (IZhou et al.l 1201 5t IZou et all 120151) 
and hence objects will have different ( u - g) o colors. In order 
to get an optimal classification for BHB/BS stars we deter¬ 
mine our boundary line using only bright stars with accurate 
photometry, namely 16 < uq < 21 and cr(u,g,r) < 0.1 mag. 
We divide the data along the (g - r) o direction, with widths of 
0.02 mag and 0.05 mag. This gives us 18 useful slices which 
we then fit using two Gaussians, 


-K2? ' 

(3) 

The fitting procedure is illustrated in the right panel of Figure 
[2] where we show bin widths of 0.05 mag. The two separate 
populations can clearly be discerned, with BHB stars lying 
to the red (i.e. right) side of the distr ibutio n. By compar¬ 
ing this figure to the similar one from IDeason et al.l (1201 ll) . 
it is evident that the SCUSS M-filter is much more gravity- 
sensitive, as the BHB/BS populations exhibit less overlap and, 
in a number of (g - r)o slices, there are clear minima in the 
(u - g )o distributions. The division is not perfect, meaning 
that BHB samples will be incomplete and contaminated by 
BSs (see Section 0, but it is clear that the SCUSS u-filter 
performs better in this regard compared to the SDSS M-filter. 

The fits for all 18 slices are shown in the left panel of Fig- 
ure|2l where the red points denote the centers of the Gaussians 
(yUi and /ii), corresponding to the centers of the BS and BHB 
populations, respectively. We adopt the median of these two 
values (i.e. H = (n i + green points in Figure0 as the 

boundary between the two populations. We have chosen p 
to represent the boundary line because this quantity is magni- 


N - N i exp 


-[(u-g)o-H 


2 \ 


+N 2 exp 


- [(u - g)o 


2cr 2 2 


8 https://www.sdss3.org/dr9/imaging/other_info.php 


tude independent. For a given (g-r) 0 slice the ratio of BHB to 
BS stars (N 1 /N 2 ) will vary as a function of magnitude, which 
means that the overall shape of the (m - g) 0 distribution will 
alter; however, even though the shape will change, it can be 
seen that // is independent of magnitude. In order to calcu¬ 
late our expression for the boundary line, we fit the values of 
/r using the following polynomial function, as shown by the 
yellow curve in Figure[2j 

(u-g\ = 1.171-0.888(g-r)„-1.531(g-r) 0 2 +11.791(g-r) 0 3 

(4) 

where stars on the right side of this curve are taken to be BHB 
stars. 

We now apply this classification to the whole magnitude 
range (16 < no < 22.5), taking all stars with errors in ugri of 
less than 0.2 mag. As mentioned above, even though Equation 
(0 was derived for stars with 16 < mq < 21, we are able to 
apply it now to the full magnitude range because /./ should 
be independent of ma gnitu de. BHB distances are calculated 
using Equation (7) of lDeason et al.l (1201 ! ). which presents a 
relation between absolute-magnitude and (g - r) a color with 
M g » 0.5 mag. 


4. BHB STARS AROUND PISCES 

Given this large and relatively clean sample of BHBs, we 
now proceed to investigate the Pisces overdensity. 

We take BHBs according to the boundary defined in Equa¬ 
tion 0, but in order to further reduce contamination we reject 
all stars within 0.02 mag of the boundary and also those which 
are more than 0.15 mag away, i.e. we only retain stars with 
0.02 < A(m - g)o = (11 - g )0 - (u - g)b < 0.15 mag. The latter 
cut does not remove many stars, but the former is important 
for reducing the number of BS contaminants in the sample. 
Note that although the error cut is 0.2 mag in each band, our 
final errors are considerably smaller. For example, the median 
error on (m - g) for our entire sample is 0.032 mag. Even at 
faint magnitudes the errors are reasonable, with BHBs at 100 
kpc having median error of 0.1 mag. 

Since our photometric precision is good and the SCUSS 
M-filter is more sensitive to gravity than the SDSS filter, the 
level of BS contamination is low. However, as can be seen 
from the right panel of Figure 0 where it is clear that the 
BHB and BS Gaussians overlap, this is not negligible and 
must be estimated. For the typical distances we wish to in¬ 
vestigate (i.e. heliocentric distances of 35 to 90 kpc, corre¬ 
sponding to 19.5 < Mo < 21.5 mag), we find that the fraction 
of BSs contaminating our sample is between 8% and 18%, 
depending on (g - r) 0. As there are no known structures in 
the foreground of the Pisces region, these BSs should be uni¬ 
formly distributed and therefore this level of contamination 
is unlikely to be problematic. Note also that our sample will 
not be complete, as a number of BHBs lie to the left of the 
boundary line. We estimate that our BHB sample is around 
80% complete. 

Now we have our sample of BHBs, we use a standard K- 
nearest-neighbor algorithm to locate the density peaks. This 
technique simply calculates the density at a point by averaging 
over a volume which includes the K nearest neighbors, where 
the choice of K depends on the problem at hand (in effect this 
is a smoothing length). For our work we have chosen a value 
of K= 10, although our results are similar if one takes K=8 or 
12. When calculating the density we have subtracted a smooth 
background model, using a flattened (q =_0.7) double power- 
law profile fit to the RR Lyrae data of I Watkins et al] (120091) ; 
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Fig. 2. — Color-color selection for BHB stars with 16 < no < 21 and <r(u,g,r) < 0.1 mag. Left panel: (u — g) o vs (g — r )o diagram for A-type stars. Red 
diamonds are the centers of BHB/BS Gaussians and green diamonds are the mid-point between these values. The yellow curve is our fit to these green diamonds. 
Stars on the right side of this yellow curve are considered to be BHB stars.Right panel: An example of the two-Gaussian fit for each slice in (g — r) o. Black 
histograms are the (if — g)o distribution for all A-type stars in each (g - r)o slice, red curves are Gaussian fits to the BHB/BS populations and the blue curves are 
the sum of these two Gaussians. 


the fit parameters are given in iFaccioli et al.l (12014 2nd row 
of Table 4). Since this profile corresponds to RR Lyrae and 
not BHBs, we have increased the normalization by a factor 
of 2, which we chose in order to match our observed distri¬ 
bution of BHBs. Given this smooth background model, we 
can calculate the significance of any volume in our 3D space, 
calculating the difference between the observed and predicted 
number of BHBs. 

We present our map of the BHB density distribution in the 
left panel of Figure [3] Since the Pisces overdensity was dis¬ 
covered in the SDSS equatorial Stripe 82 data, we first in¬ 
vestigate that region (i.e. | Dec | < 1.25°). Indeed there 
is a strong overdensity located around -5° < RA < -15°, 
which is precisely where it was previously identified. The 
other features of note in this region are the additional struc¬ 
tures at -25° < RA < -20°, but these are most-lik ely the 
edge of the Hercules-Aquila Cloud dBelokurov et alJ 120071 : 
Simion et al.l2014l). which was also detected in Stripe 82 (e.g. 
Watkins et al .1120091) . 

When we compare the distance to our Pisces detection 
(middle panel of Figure [3j again restricting ourselves to just 
the Stripe 82 region), we find a weak secondary detection 
behind the main clump at around 95-100 kpc. This more- 
distant feature has only 5 members and hence may be a spu¬ 
rious detection, but the significance is noticeable (at around 
2cr) because the model predicts very few stars this far out in 
the halo. The detection is unlikely to be caused by extra- 
galactic contamination; XDQSO classifies all 5 members as 


stars and so the overdensity is unlikely to be caused by back¬ 
ground quasar contamination and, since r ~ 20.5 mag for 
these stars, the SDSS star-galaxy classification should be ro¬ 
bust dLupton et al.l200ll) . This detection cannot be a ‘shadow’ 
of mis-classified BSs from the main structure, since the offset 
in distance modulus is around 0.6 mag, not the 2 mag that one 
would expect if these were misclassified BSs. Also, with the 
exception of one star, all are far from the BS/BHB boundary 
and so this reinforces our belief that these are not BS contam¬ 
inants. 

If one compar es o ur distance distribution to the RR Lyrae 
detection from I Watkins et all (120091 Figure 16), then one also 
sees a hint of bimodality and a similar spread in distances 
(around 40 kpc). It should be noted that the peak of our BHB 
detection is slightly offset from the RR Lyrae detection, with 
distances of 75 and 80 kpc, respectively. Since there is more 
uncertainty in the BHB distance calibration (as opposed to 
the RR Lyrae distances), it may be that this is out by 0.1 to 
0.2 mag. This could be due to deficien cies in our adopted 
absolute-magnitude relation ( Equatio n 7lDeason et al.li20lTl) . 
such as a metallicity bias (iFermani et al.ll2013l) . or differences 
in the adopted extinction values. 

The next step is to expand our search beyond the Stripe 82 
region. As can be seen the left panel of Figure [3j the wide 
coverage of the SCUSS data, combined with the clean sepa¬ 
ration of BHBs from BSs, allows us to tentatively detect an 
extension of the overdensity to both lower and higher de¬ 
clinations. These manifest themselves as separate clumps. 
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which seem to be connected to the peak overdensity lying 
at RA = -10° and Dec = -3°. The significance of these 
two regions is 3.5 and 3.0cr, respectively, for the lower- and 
higher-declination regions. There are some additional clumps 
at -25° < RA < -20°, but again these are most-likely the 
edge of the Hercules-Aquila Cloud. 

It seems that the Pisces overdensity is extended like a 
stream, passing from the bottom-left to top-right of the fig¬ 
ure. In order to confirm this we now analyze the distances of 
the BHBs in this feature. We define a path for the stream (de¬ 
noted Ap and Bp; see Appendix) and take BHBs within 2 deg. 
The density map of these is shown in the right panel of Figure 
[3] Here we can see that the detections are aligned beautifully 
along a distance gradient, confirming our interpretation of a 
stream. Again the curious cloud of distant BHBs is noticeable 
beyond 100 kpc, but is clearly unrelated to the stream. The 
distance gradient along the stream is around 1.3 kpc per deg, 

D P = 1.3 X A P + 72, (5) 

where Dp is the helio-centric distance of the stream in kpc. 

We calculate the significance by taking a volume encom¬ 
passing our Pisces detection (-12 < Ap < 10 deg, |Bp| < 
2 deg, |D - Dp| < 10 kpc) and calculating the number of 
BHB stars predicted from the smooth model, then compare 
this to the observed number. The model predicts that there 
should be 8-9 stars in this volume, while we actually find 27 
stars. This implies that our detection is 6.2cr. The signifi¬ 
cance varies depending on how one defines the background 
model, for example the flattening or the normalization, but 
the variation is small and so we conclude that the signif¬ 
icance is around 5-lcr. The significance of the new de¬ 
tections that extend away from the central region, i.e. at 
(-12 < Ap < -5 deg) and (3 < Ap < 10 deg), is 3.5 and 3.3cr, 
respectively. 

Although we find that 27 stars belong to our detection, this 
is only an approximate number. As discussed above, some of 
these may be BS contaminants, while some bona fide BHBs 
may lie beyond our BHB/BS boundary line. An additional 
four objects were rejected by our quasar cut (Equation^, of 
which two are likely to be stars according to XDQSO. On the 
other hand, of the 27 objects that passed our quasar cut, one 
was classified as a quasar by XDQSO. 

There are a couple of notes of caution which should be ad¬ 
dressed. Firstly, there appears to be a gap in the stream at 
Kp ~ 4 deg and Dp ~ 65 kpc. We have checked whether this 
could be due to bad fields or patchy extinction, but neither are 
found at this location. However, the density of BHBs is low 
and so this gap could just be due to statistical fluctuations. 
The other cautionary point is the clump of material around 
-14° < RA < -8° and -10° < Dec < -5°. Although this 
is statistically significant (at around 3cr), the distance distri¬ 
bution is very broad and there is no obvious clumping at the 
distance of Pisces. 

We have chosen to focus on the region around the Pisces 
overdensity, even though the SCUSS footprint is much larger. 
However, if we extend the area to the entire footprint, we do 
not detect any significant new structures. 

5. GIANT BRANCH EXPLORATION 

In an attempt to compliment this analysis, we have also an¬ 
alyzed data obtained from the 3.6m Canada-France-Hawaii 
Telescope (CFHT) using the 1-square-degree field-of-view 
MegaPrime/MegaCam camera. We chose 8 min exposures 


in r (seeing of ~0.5") and 12 min in g (seeing of —0.9"), al¬ 
lowing us to reach around 24 mag in both bands. Since we 
could not contiguously map the full region around Pisces, we 
limited ourselves to 30 pointings (i.e. 30 sq. deg.) spread 
over a wide area. The location of these fields are shown in the 
right panel of Figure [4] Although the data reach to 24 mag, 
the star-galaxy separation (which is dependent on the seeing 
and signal-to-noise ratio) becomes problematic below around 
r « 22 mag. In order to alleviate this problem we cross-match 
with SCUSS u band and apply the following mask to isolate 
the stellar locus (see, for example [Strateva et al.ll200lh 

| (u - r) o - 0.37(g - r) 0 | < 0.75, (6) 

and also apply (u - r) o > 0.75 to remove quasars. 

The resulting Hess diagram for our data is shown in the left 
panel of Figure [4] The main-sequence turn-off for the stellar 
halo is clearly visible, but there is no obvious detection of any 
sequence belonging to the Pisces overdensity. To guide the 
eye, we have included in this figure a box corresponding to 
the isochrone for an 8 Gyr population with [Fe/H] — —1.5 
dex, located at 75 kpc. To give the box a finite size we have 
shifted it by +5 kpc and by +0.03 mag in (g - r) q. Due to the 
limitations of the SCUSS data, completeness begins to drop 
around r ~ 22.5 mag and we are therefore unable to detect 
the main-sequence turn-off. If the giant branch is present, it 
is not immediately obvious although it may be hidden behind 
the main-sequence of the stellar halo, which is considerably 
more dense. 

We investigate further by calculating the number of stars 
within the isochrone box shown in Figure [4] In order to ac¬ 
count for variations in the background density of stars across 
our fields, we need to normalize this number. We do this by 
dividing the number of stars inside the isochrone box N m by 
Nout, where N out is the total number of stars in the same mag¬ 
nitude range (22.5 < tq < 20.5, 0 < (g - r)o < 0.8) ex¬ 
cluding stars inside the isochrone box and those within 0.1 
mag in (g - r) o of the box boundary. We calculate this frac¬ 
tion Nj n /N out as a function of Bp, the cross-stream angle in¬ 
troduced in the previous section and plot this in the middle 
panel of Figure 0 At most locations this fraction is relatively 
constant at around 0.425, but for two bins close to the centre 
of the stream the value is notably above this (at a significance 
of 1- to 2-cr, where errors are assumed to be Poissonian). 

Finally we split our data into the individual CFHT fields. 
We first estimate the background by averaging all Bp bins ex¬ 
cept the two at Bp = -2° and 0°. We then calculate this frac¬ 
tion for each of our fields and, in the right panel of Figure [4] 
show with filled boxes the fields for which the fraction is at 
least 0.5cr above the background value. From this figure one 
can see that there is a tendency for fields with larger fractions 
to lie within 2 deg. of the proposed stream plane (shown by 
the dashed lines). However, this is inconclusive, with most 
fields having excess at only 1- to 2-cr. 

In summary, our search for the Pisces giant branch shows 
tentative evidence for alignment with the stream plane in the 
previous section, but the results are inconclusive and veloci¬ 
ties would be required to clarify the situation. 

6. DISCUSSION 

In this paper we have provided an extended view of the 
Pisces overdensity, utilizing deep u band data from the 
SCUSS survey. With the SCUSS u band data, which goes 
around 1 to 1.5 mag deeper than single-epoch SDSS, we have 
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RA (degree) RA (degree) A P (degree) 


Fig. 3.— The density of BHBs around the Pisces overdensity, calculated using an 10-nearest-neighbor method. The original detection was in the narrow Stripe 
82 region (dotted lines in the left panel) at —5° < RA < —15° and the corresponding distances to our BHBs overdensities in this region are shown in the middle 
panel. The left panel shows the BHB density for areas surrounding the original detection. The overdensity shows an extension along a stream, denoted by the 
black solid line (see Appendix for the coordinate transformation). The stream plane is defined by two angles (Ap and Bp), where Ap is oriented along the stream 
and increases with decreasing RA, and Bp is oriented across the stream and increases with increasing Dec. The dashed lines lie 2 deg either side of this path 
and the distances of BHBs within this region are shown in the right panel. The fact that the distances lie along a well-defined distance gradient confirm our 
interpretation that Pisces is part of a stream. 


0.0 0.2 0.4 0.6 

( a - r )o 





Fig. 4.— Our search for the giant branch of the Pisces overdensity, using data from CFHT. The left panel shows a Hess diagram for all of our fields, with the box 
corresponding to an isochrone for an 8 Gyr population with [Fe/H] = —1.5 dex, located at 75 kpc. To give the box a finite size we have shifted it by ±5 kpc and 
by ±0.03 mag in (g - r) o. The middle panel shows the number of stars within the isochrone box, normalized by the total number in the same magnitude range 
(22.5 < ro < 20.5, 0 < (g — r)o < 0.8) but excluding stars inside the isochrone box and those within 0.1 mag in (g — r )o of the box boundary. This is shown as a 
function of cross-stream angle Bp and errors are assumed to be Poissonian. Finally, in the right panel we show the location of our CFHT fields. Filled squares 
are those which have fractions greater than the background level by at least 0.5<r. 


used BHB stars as tracers and analyzed the distribution of 
overdensities for ~1000 sq. deg. around the original Pisces 
detection. We have found that Pisces appears to be part of a 
stream, with a clear distance gradient (Figure [3]). This stream 
is around 5 deg. wide and we appear to trace it for 25 deg. in 
length. Given the considerable distance to this stream (60 to 
80 kpc) this makes it the largest structure (by volume) in the 
outer halo after the tidal tails of the Sagittarius dwarf galaxy. 
We find that there are 27 BHBs associated to our extended 
detection of the overdensity, although we cannot report this 
number with any high degree of accuracy due to three effects: 
firstly, some of these BHBs will not belong to the overden¬ 
sity and are either smooth halo BHBs or foreground BSs; sec¬ 
ondly some BHBs will be lost due to incompleteness, espe¬ 
cially those which have smaller (bluer) (u - g)o and are hence 
indistinguishable from BSs; and thirdly, there may be mem¬ 
ber stars outside the Ap range we have detected. We can 
roughly estimate the influence of the first two effects, using 
the contamination and completeness fractions calculated in 
Section [4] After accounting for these, we predict that total 
number of BHBs belonging to our Pisecs detection is around 


21 stars. With this rough number in hand, we can speculate 
as to the progenitor of the stream. Assuming we have de¬ 
tected the entire stream and that the progenitor has fully dis¬ 
rupted, then this number of BHBs suggests that the progen¬ 
itor was likely somewhere between a smaller classical dwarf 
spheroidal such as Draco or Sextans (e.g. lAparicio et al.l2001t 
iLee et alJl2003l). and a larger ultra -faint dwarf like Canes Ve- 
natici I (e.g. lOkamoto et al.112012i) . 

Compared to the other Pisces detections^ which are 
restricted to the Stripe 82 region dSe sar et al.1 l2007t 

iKollmeier et al.ll2009t IWatkins et al1l2009USesar et*akll2010l) . 

our detection of a tangential extension is new. Before, the 
only ten tativ e evidence of an extension was reported by 
ISharma et al, ' d2010l) . Using M-giant tracers and a group 
finding method, they report that the Pisces overdensity ex¬ 
tends from -20° <RA< 25°, 1.25° <Dec< 25°, with a dis¬ 
tance of 103+51 kpc (see Figure 9 of their paper). Com¬ 
pared to our detection, there is a similarity in the region of 
-20° <RA< -10°, 1.25° <Dec< 11°, but their detection 
at RA> -10° is totally different from ours; here we find 
the overdensity extending to negative declinations, but the 
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ISharma et al.l detection extends to larger positive declinations. 
Furthermore, this detection does not overlap with the Stripe 
82 detection. A mo re-recent M-giant map of that region is 
presented (iDeason et al.l 12014 ). Although not remarked on 
by the authors themselves, an inspection of the left panel of 
their Figure 5 shows an overdensity at the same location as 
the Pisces detection, seemingly oriented in the same direction 
as our BHB detection. 

The only discrepancy between our results and the exist¬ 
ing Stripe 82 detections is the offset in distance, with ours 
being around 10 % smaller than those estimated from RR 
Lyrae. We ar e una ble to determine the cause of this offset. 
Like lWatkins et all (120091) we also detect a hint that the Stripe 
82 detection may be split into a near and far component, but 
further work needs to be carried out to confirm this. In par¬ 
ticular is it unclear whether this could explain the multiple 
kinematic groups in t he S tripe 82 detection , as claimed by 
iKollmeier et al.l (120091) and lSesar et al.l (120101) . It appears that 
the distanc es to bo th fe ature s are similar (see, for example, 
Figure 1 of lSesar et al.ll2010l) . but with such small numbers of 
stars it is hard to make firm conclusions. It is also clear that 
our detection of the stream is rather ‘lumpy’, which may re¬ 
flect the lumpiness of the stream itself, or possibly due to the 
sparseness of the tracer population. 

Since we have detected a stream, it is natural to ask whether 
the core of the progenitor remains intact. There is a dwarf 
galaxy in the vicinity of the stream, the recently detected 
Pisces II dwarf (RA=344.6°, Dec=5.9°) ([Belokurov et al.l 
l201d ). If we transform its sky coordinates to our stream coor¬ 
dinates, we find that it has Ap = 7.9° and Bp = 0.7°. However, 
despite it’s close proximity (in projection) to the stream plane, 
its distance of 180 kpc precludes it from being the progenitor, 
at least not if it is part of the same wrap of the stream. Given 
our detected distance gradient, it is now possible to model the 
orbit of the stream and embark on a large scale hunt for the 
progenitor. 
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APPENDIX 

THE COORDINATE TRANSFORMATION 

Here, we provide the equations for converting from Equatorial (a, d) to the Pisces Overdensity coordinate system (Ap, Bp). 
The orbital pole of the Pisces overdensity is set to (ap, dp), and its center is set to (a c , 6 C ). Details of the derivation are shown 
below. 

Spherical coordinates can be converted to a right-handed Cartesian coordinate system using, 


X 

f cos a cos 6 ' 

Y 

= sin a cos 6 

u 

V sin d 


The orbital pole of the Pisces system (ap, dp) can be rotated into the new Cartesian coordinate system as follows. 


' X' ’ 


( R, x 

Y' 

= 

Ry Y 

U' 


, R- z , 
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where 


With X',Y' and Z , 


A i xA 2 
IAi xA 2 r 




[cos ap cos dp, sin ap cos dp, sindp], A 2 — [0,0,1] , 
R x x Ai 

Ry = | | ’ 


i?, = [cos ap cos dp, sin ap cos dp, sindp]. 


A P = atan2(Y,X) 


and 

Bp = arcsin(Z ) , 

where tan(atan2((Y',X')) = (Y'/X'). 

If the center of the Pisces overdensity is set to ( a c , 6 C ), then we need to put this in phase with the original coordinate system, 

Ap = Ap - Ape 


Bp = Bp - Bpc 


where Ape and Bpc is for (a, <5)= (a c , 6 r ). 

Our orbital plane uses (ap, 6p) = (19°, 47°) (a c , 6 C ) = (-10°, 0°), which results in the following transformation, 

Ap = atan2(-0.13954893 cos a cos 6 - 0.71791667 sin a cos 6 + 0.68199837 sinb, 

0.98162711 cos a cos 6 - 0.19080906 sin a cos d) + 0.73139161° (Al) 

Bp = arcsin(0.13013147 cos a-cos <5 + 0.66946810 sin a cos d + 0.73135370 sin d) 


Note that in this system Ap is oriented along the stream and increases with decreasing a, and Bp is oriented across the stream 
and increases with increasing d. 

The reverse transformation from the Pisces Overdensity coordinate system (Ap, Bp) to the Equatorial system (a, d) is, 

a = atan2(-0.19080907 cosAp cos B P -0.71791661 sin Ap cos B P + 0.66946810 sinB P , 

0.98162729 cos Al, cos B P - 0.13954890 sin Al, cos B P + 0.13013147 sin B P ) 

, P (A2) 

6 = arcsin(0.68199837sinApCosB P + 0.73135370 sinB P ) 

Ap = A P -0.73139161° 
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